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^ Abstract: A low-background germanium spectrometer has been installed and is being operated 

in an ultra-low background shield (the Gator facility) at the Gran Sasso underground laboratory 
in Italy (LNGS). With an integrated rate of ~0. 16 events/min in the energy range between 100- 
2700 keV, the background is comparable to those of the world's most sensitive germanium detec- 
ts tors. After a detailed description of the facility, its background sources as well as the calibration 
and efficiency measurements are introduced. Two independent analysis methods are described 
and compared using examples from selected sample measurements. The Gator facility is used to 
^ screen materials for XENON, GERDA, and in the context of next-generation astroparticle physics 
O facilities such as DARWIN. 
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1. Introduction 

Gamma-ray spectroscopy offers a standard method for material screening and selection for rare- 
event seaixhes, such as the direct detection of Weakly Interacting Massive Particles (WIMPs) or 
the search for the neutrinoless double beta decay [0]. An ultra-low level germanium spectrometer in 
a dedicated low -background shield (the Gator facility) has been built and installed at the Laboratori 
Nazionali del Gran Sasso (LNGS), Italy. While the facility is being operated mainly in the context 
of the XENON program [0] [§] [0] [0], recendy it also has been used for the GERDA experiment [H 
as well as for R&D purposes within low-background particle astrophysics. 

This paper is organized as follows: in the next section, the Gator facility is described in detail. 
In the third section, calibration measurements, efficiency determinations for the different sample 
geometries, and cross-checks with standard, calibrated samples are described. In the fourth section, 
two different data analysis methods are introduced. The fifth section describes the determination 
of the main background sources of the facility as well as selected results for sample measurements. 
In the final section, a summary and discussion of future plans are presented. 
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2, The Gator facility 



The design of the faciUty has been inspired by the layout of the world's most sensitive germanium 
spectrometers, operated at LNGS in connection with the Borexino and GERDA experiments [1I][H]. 
The core consists of a high-purity, p-type coaxial germanium (HPGe) detector with 2.2 kg of sen- 
sitive mass with a relative efficiency of 100.5% ^ The detector construction has been performed in 
close cooperation with CanbeiTa semiconductors [ITTll . using only materials with ultra- low intrin- 
sic radioactive contamination. The cryostat is made of ultra-low activity, oxygen-free copper with 
the cooling provided by a copper dipstick ('cold finger') in thermal contact with a liquid nitrogen 
bath. The cryostat is of U-type, with the cooling rod shaped in a right angle below the cryostat to 
avoid direct line-of-sight to the outside (see Figure [I]). While the field-effect transistor (FET) used 
for the charge readout is cooled and close to the detector, the preamplifier is placed outside the 
low-background shield, since it contains more radioactive components. 

The shield of the detector has been designed to provide a large sample capacity, an ultra-low 
background and easy access to the germanium spectrometer itself. The sample chamber, with a 
dimension of 25x25x33cm^, is surrounded by 5 cm (7 cm for the base plate) of oxygen- free, 
radio-pure copper from Norddeutsche Affinerie WTSl . Residual surface contaminations of the cop- 
per plates were removed by treating them with diluted sulfuric and citric acid solutions, followed 
by cleaning with deionized water. All steps were performed under clean-room conditions. The 
copper is surrounded by four layers of lead from Plombum [lOll . each 5 cm thick. The inner 5 cm 
lead layer has a nominal ^'°Pb activity of 3 Bq/kg, while the outer 3 layers were built from lead 
with a higher ^^°Pb activity of 75 Bq/kg. All lead bricks were cleaned with ethanol before being 
installed in the shield. Their aiTangement is such that no direct line-of-sight to the HPGe detector 
is possible. Lai^ge copper plates, which close the sample cavity and caiTy the upper lead bricks, 
were placed on sliding rails allowing for easy opening and closing of the shield: The full sample 
chamber can be accessed when the top shield is open. The lead is surrounded by 5 cm of polyethy- 
lene as a shield against ambient neutrons, and the entire shield is enclosed in an airtight aluminum 
box in order to prevent radon entering the system. 

Figure |I] shows a schematic view of the detector and shield configuration. A sample handling 
and glove box made of plexiglass, including an airlock system, is placed on top of the aluminum 
housing. The entire system is continuously flushed at slight over-pressure with boil-off nitrogen 
gas to suppress radon diffusion into the shield. A 5 mm inner diameter PTFE tube allows sealed 
calibration sources to be brought close to the germanium detector. The samples to be screened are 
first cleaned in an ultra-sonic bath of ethanol (where applicable), then enclosed in a sealed plastic 
bag for transportation to the underground site, after which they are stored for a few days under 
nitrogen atmosphere above the closed chamber. This allows trapped radon and plate-out ^^'^Rn and 
^^^Rn progenies to decay before the actual counting starts. 

To remotely check the stability of the detector with time, a monitor system has been installed. 
The liquid nitrogen level, the flow of nitrogen gas, the leakage current across the germanium diode 
and the high- voltage applied to the diode are read every 2 minutes, while the overall trigger rate is 
read every 6 hours. The liquid nitrogen level is measured with a level-meter consisting of two 40 cm 

' The quoted efficiency is defined relative to a 7.62 cm x 7.62 cm Nal(Tl) crystal, for the 1.33 MeV *°Co photo-peak, 
at a source-detector distance of 25 cm IIIUII . 
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Figure 1. Schematic figure of the Gator facility at LNGS. The HPGe detector (a) with its cold finger (b) and 
dewar (c), the open sample chamber, the copper (d) and lead (e) shield with the sliding door, the glove box 
(f) and polyethylene shield (g) can be seen. 

long concentric metal tubes acting as a capacitor whose capacitance is read out with a universal 
transducer interface board. The nitrogen gas flow into the shield is monitored and regulated with 
an electronic flowmeter and the leakage current is measured as the voltage drop across the feed- 
back resistor and read out with an analog-to-digital converter. The relevant parameters are plotted 
versus time and can be monitored on a web interface which is refreshed every 10 minutes. In case 
pre-defined thresholds for these parameters are not met, email and SMS alarms are being issued. 
The data acquisition consists of a high- voltage unit, a spectroscopy amplifier, and a a self triggering, 
16 k channel multichannel analyzer. The spectra, acquired and cleared every 6 hours, are saved in 
ascii files and analyzed offline. 

Prior to its installation at LNGS, the germanium spectrometer was operated in the Soudan 
underground laboratory in northern Minnesota, within the SOLO facility W\M . At the Soudan labo- 
ratory, the detector was used for the screening of XENON 10 ED materials, and several background 
runs were acquired [1T511 . During its water- and ground-based transportation to LNGS, the detector 
was exposed for several months to the cosmic ray flux at the Earth's surface, leading to cosmic 
activation of the crystal and the surrounding copper of the cryostat. The shield was improved in 
October 2008, and it was cleaned once again in February 2009. 

Table |I] shows the integral background in the 100-2700 keV region in the Soudan configuration, 
and for three measurements taken at LNGS. The integral counting rate at LNGS has been constantly 
decreasing due the decay of cosmogenic radio-nuclides such as ^^Mn, ^^Co, ^^Co and ^^Zn, with 
typical half-lives around one year on the one hand, and due to an improved shield and overall 
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sealing of the system on the other hand. This is reflected in the decrease of prominent lines from 
^'^Bi and ^^"^Pb, as shown in Table It gives the integral counting rates in the ib3a-regions for the 
main primordial 7-lines, as well as for the '^^Cs, ^"Co and lines, along with a comparison with 
the GeMPI detector, which is one of the world's most sensitive low -background spectrometers [ED. 



Table 1. Integral background counting rates for Gator as measured at Soudan and at LNGS in three different 
runs. The integral is evaluated in the energy range [100, 2700] keV. 



Run 


Lifetime [days] 


Rate [events/min] 


Gator at Soudan 


22.96 


0.842 ± 0.005 


Gator at LNGS (09-2007) 


14.90 


0.258 ± 0.003 


Gator at LNGS (10-2008) 


22.59 


0.186 ±0.003 


Gator at LNGS (04-2010) 


51.43 


0.157 ±0.001 



Table 2. Background counting rates (in events/day) in the ±3(7-regions for the main primordial and the 
gamma lines of '^^Cs, ^"Co and "^^K. 



Energy 


Chain/nuclide 




Peak integral background rate [counts/day] 




[keV] 




Gator 


Gator 


Gator 


Gator 


GelVIPI \M 






(Soudan) 


(LNGS, 
09-2007) 


(LNGS, 
10-2008) 


(LNGS, 
04-2010) 




239 


232Th/212pb 


1.1 ±0.7 


0.7 ±0.1 


0.13 ±0.08 


<0.5 


NA 


911 


232Th/228Ac 


0.9 ± 0.3 


0.4 ± 0.2 


0.4 ±0.1 


<0.5 


<0.2 


352 


238u/214pb 


4.9 ± 0.7 


4.3 ± 0.7 


1.1 ±0.2 


0.7 ± 0.3 


<0.5 


609 


238u/214Bi 


4.5 ± 0.5 


4.0 ± 0.5 


1.1 ±0.2 


0.6 ± 0.2 


0.50±0.45 


1120 


238u/214Bi 


1.6 ±0.3 


2.7 ± 0.4 


1.3 ±0.2 


0.3 ±0.1 


NA 


1765 


238u/214Bi 


1.3 ±0.2 


1.5 ±0.3 


0.2 ±0.1 


0.08 ± 0.06 


NA 


662 


137CS 


2.9 ± 0.4 


0.5 ± 0.3 


<0.5 


0.3 ±0.1 


NA 


1173 


6OC0 


0.5 ±0.1 


0.5 ± 0.2 


0.5 ± 0.2 


0.5 ±0.1 


0.6±0.4 


1332 


6OC0 


0.6 ±0.1 


0.6 ± 0.2 


0.5 ±0.1 


0.5 ±0.1 


0.4±0.3 


1461 


40k 


5.8 ± 0.4 


0.5 ± 0.2 


0.4 ±0.1 


0.5 ±0.1 


0.6±0.4 


2615 


208 rp J 


0.7 ±0.1 


0.2 ±0.1 


0.2 ±0.1 


0.2 ±0.1 


NA 



Figure (top) shows the comparison between the latest background spectrum acquired at 
LNGS (2010), a spectrum taken in the SOLO facility at Soudan (2007) and the background of the 
GelVIPI detector It also shows (bottom) the Gator spectrum underground at LNGS prior to its 
installation in the shield, inside the shield, and inside the shield with the radon protection system 
on. The background decrease is more than four orders of magnitude, and shows that a careful 
shielding is needed even when the detector is operated deep underground. 
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Figure 2. (Top) Background spectra of Gator at Soudan (red), at LNGS (black) and the spectrum of 
the GeMPI detector [§] (blue). (Bottom) Gator background spectrum at LNGS: outside the shield (black), 
inside the shield (blue) and inside the shield with the radon protection system on (red), clearly showing the 
suppression of the main gamma lines associated with radon decays. 



3. Calibration measurements and efficiency determination 

The HPGe detector is calibrated regularly with radioactive sources such as '"^Cd, ^^^Ba, '^^Cs, 
^°Co, ^'Co, ^^Na, ^"^Mn and ^^^Th. In Figure H (left) the comparison of the spectrum obtained from 
a ^"Co calibration with the one from a Monte Carlo simulation of the source-detector geometry 
is shown. The FWHM of the two ''"Co lines at 1173 keV and 1332keV are 2.5 keV and 3.0keV, 
respectively. The energy resolution of the detector, defined here as the ratio of the a to the mean 
energy of the gamma line, is shown in Figure ^ (right). 

The calculation of the detection efficiencies of the various gamma lines used in the analysis of 
the experimental data (see Section ^ is based on Monte Carlo simulations using Geant4 [ITHl . For 
each measured sample, a detailed geometry is included into a Geant4 model of the facility. The 
efficiency £ of a specific y-line is defined as the ratio between the number of events detected in 
the line to the number of gammas of that energy emitted by the source. In order to simulate each 
decay chain, the G4RadioactiveDecay class, which takes into account the branching ratios for the 
different gamma lines in one decay, is used. 
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Figure 3. (Left) Comparison of the Gator spectrum from a "Co calibration (black, data points) with the 
one from a Monte Carlo simulation of the source-detector geometry (red, solid). (Right) Energy resolution 
(defined here as the ratio of ff/energy) as a function of energy. The solid curves represent a fit using the 
function a^{E) = E^{2.35 x 10"^) +£(7.70 x IQ-"*) + (4.43 x IQ-^). 

To cross-check the efficiency determination, a measurement of two extended sources and a 
comparison with the certified values for their activities has been performed. The sources used 
for this measurements, which had similar dimensions and weights, are CANMET-STSD2 (from 
the Canada Centre for Mineral and Energy Technology) and IAEA-SoiI6 (from the International 
Atomic Energy Agency). Both sources are soils from different places on Earth, which have been 
thermally treated and sieved several times in order to destroy any remaining organic matter. The 
homogeneity of the material is certified by the provider. Tables shows the results Gator's screen- 
ing of these two sources and a comparison with the certified values. Within the uncertainties of 
the measurements, we find a good agreement. We also show in Figure the comparison of the 
efficiencies as determined by the Monte Carlo method, with those from the data, as a function of 
energy, as well as the relative difference among these. 

These results indicate that the measurements performed with our spectrometer provide a re- 
liable value for the activity of a given sample. Further cross-check were done by using well- 
calibrated, commercially available point sources. 



4. Data analysis methods 

Two methods are used to determine the concentrations of radioactive nuclides in a given sample. 
In the first method, the most prominent y-lines are analyzed, using efficiencies as determined by 
a full Monte Carlo simulation. In the case of ^^^U, the /-lines from the daughters of ^^^Ra (^^^Pb 
and ^^^Bi) and in case of ^^^Th the /-lines from ^^^Ac and from ^^^Pb, ^'^Bi and ^^^Tl are used. In 
the second method, the overall data spectrum is compared to the one obtained from a Monte Carlo 
simulation, after subtracting the measured background spectrum, and the activities are determined 
from the best fit. The results from the two methods agree within the statistical errors, as shown in 
Section^ Both methods are explained in more detail below. 
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Table 3. Results from screening the CANMET-STSD2 (493 g) and IAEA Soil6 (530 g) sources and com- 
parison with certified values provided by two agencies (for details, see text). 



Nuclide 


STSD2 Activity [Bq/kg] 
Gator results Certified values 


228xh 
40k 


75±4 70±5 
230±30 230±10 
590±10 540±20 


Nuclide 


Soil6 Activity [Bq/kg] 
Gator results Certified values 


226Ra 


88±5 80±7 
57±2 54±2 
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Figure 4. Upper panel: Efficiencies as a function of energy for the two certified sources, as determined from 
data (filled symbols) and Monte Carlo (open symbols). Lower panel: relative difference between simulated 
(MC) and measured (D) efficiencies. 



4.1 Analysis of gamma-lines 

The first method is based on counting events at the location of the most prominent lines, after 
subtracting the background spectrum closest in time. The Compton background, estimated from 
the regions left and right of a peak, is subtracted as well. The decision on whether the signal 
exceeds the background is based on comparing the net signal number of counts 

Snet=S-B-ts/tB-Bc (4.1) 

with the so-called detection limit (the level of a true net signal that can be detected with a given 
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probability) as defined in flU for a ~95% C.L.: 



= 2.86 + 4.78 WBc + S- - + 1-36. (4.2) 



S is the number of counts in the ±3(7 -region around a peak, B and Be are the number of background 
and Compton-background counts in the same region, and ts, ts are the measuring times for signal 
and background, respectively. For each peak, three cases are considered HTEII : 

1. Snet < 0: the upper limit is set to (no net contribution from a signal) 

2. 0< Snet < L^i- the upper limit is set to S^et (there is an indication of a signal, but it can 
not be confirmed for the existing background level and sample exposure) 

3. Snet > Ld- the detection limit is exceeded (clear indication for a signal at 95% C.L.) 
For the third case, the specific activity and its 1 a error is calculated as 

A[Bq/kg] = ^^ with ^ = ^. (4.3) 
r-e-m-t A Snet 

with the peak efficiencies e as determined by Monte Carlo simulations, the branching ratio r for the 
specific line, the mass m of the sample (in kg), and the measuring time t (in seconds). For the case 
in which an upper limit is reported, S^t is replaced by Lj or by {S„et +Ld) in equation (|4.3D. As 
concrete examples, we show the above quantities in Table ^, together with the determined specific 
activities or upper limits for a copper and a stainless steel sample, using different gamma lines from 
their measured spectrum. 

Table 4. Examples of upper limit or specific activities calculation. Details are given in the text. 



Sample 


Used line [keV] 


Snet 


B-ts/tB 


Be 


Ld 


Condition 


Activity [mBq/kg] 


Copper 


239 








93 


49 


Snet ^ i-'d 


<0.33 




352 


-5 


19 


71 


48 


Snet ^ i-'d 


<0.36 




1173 


42 


6 


19 


27 


Snet ^ i-'d 


0.24±0.06 


Stainless steel 


352 


66 


7 


58 


42 


Snet > i-'d 


4.3±0.9 




1173 


236 


2 


19 


25 


Snet ^ i-'d 


7.2±0.9 




1461 


-3 


3 


10 


21 


Snet ^ i-'d 


<5.7 



4.2 Fit of the data to a simulated spectrum 

The individual activities are also determined by a global fit based on a ;^^-minimization of the 
simulated spectrum to the experimental data. The aim is to model y, the measured number of counts 
in each energy bin, as a function of x, the energy bin value from the Monte Carlo simulation, with 
a functional dependence of the form 

M 

y=Y.ak-fk{x)^ (4.4) 

A:=0 
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where M is the number of simulated radioactive isotopes, >0 are the scaling factors for each 
isotope, which are kept as free parameters, and fk{x) are M Monte Carlo spectra which are already 
smeared with the energy resolution of the detector. For uncorrected statistical errors in each energy 
bin /, the X'^ value is defined as 

to 

where (T; = is the variance in the observed number of counts in each bin and N is the number 
of energy bins over which the fit is performed. The statistical uncertainty in the Monte Carlo 
component is negligible. Before the analysis, the measured background spectrum is subtracted 
from the sample spectrum. The outcome of this procedure are the scaling factors it for 
every decay chain/isotope k. In the subsequent analysis, these factors are assumed to describe a 
normalized Gaussian probability distribution function (PDF) with mean qq = au and a = Aa^;. 

The method applied to decide between a real detection and an upper limit is based on fUTIl . 
If the lower of the two symmetric limits providing 95% statistical coverage is positive, a peak 
detection is claimed and the activity is calculated from Uk, taking into account the measuring time 
t in seconds, the sample mass m in kg, and the number of simulated events risim- The relative it la 
eiTor is given by Aayt/a^.: 

Bq/kg = — with — = . (4.6) 

m t Ak Uk 

If the lower hmit is equal or less than zero, no detection can be claimed at 95% C.L. and an upper 
limit is given. Its value is determined by the 95% quantile of the positive part of the Gaussian 
PDF defined by = and a = Aajt. The limit on the activity (95% C.L.) is calculated using a^p 
in equation ([4.61). 

As an example. Figure shows a measurement of a stainless steel sample: the data spectrum 
is compared with the Monte Carlo sum spectrum, and the individual contributions from ^^^U, ^^^Th, 
"^^K, and ^"Co, as given by the best-fit, ai^e shown. The derived activities and upper limits for this 
sample ai^e given in Table ^ 



5. Results 



In this section, we first outline results obtained from a detailed study of the background of the 
facility, using the analysis method described in Section We then present screening results 
for a few selected samples, the goal being to compare the outcome of the two analysis methods 
introduced in the previous section. Results from a much larger selection of screened samples are 
given and discussed in detail in [EDU . 

5.1 Background Analysis 

To model the residual background of the Gator facility, the detailed geometry of the crystal, cryo- 
stat system and shields has been simulated with Geant4. The following potential contributions to 
the background have been simulated: the natural decay chains of ^^^U, ^^^Th and "^^K decays in 
the copper of the shield and of the cryostat, the decays of the cosmogenic radio-nuclides ^^Mn, 
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Figure 5. Best-fit of the Monte Carlo simulations to the measured data for a stainless steel sample. Data 
is shown with error bars (red), the full spectrum as determined by the simulation is shown as the grey solid 
curve. The individual contributions from the decays of ^Co (blue), ^^**U (green), ^^^Th (magenta), ^^Mn 
(cyan) and ^^Co (brown) are also shown, along with the measured background spectrum (black). 



^^Zn,^^Co, ^^Co, ^"Co in the Ge crystal, the Cu of the cryostat and of the shield, ^^°Pb decays in 
the innermost Pb layer, and ^^^Rn decays inside the shield.^ Other materials in the cryostat (PTFE, 
mylar and kapton) were neglected in the simulation due to their much lower mass compared to the 
copper (a few tens of grams and milligrams, respectively). 

Figure ^ (top) shows the comparison of the latest measured Gator background to the sum spec- 
trum obtained from the Monte Carlo simulations. On the bottom, the individual contributions of 
the main simulated background components are shown. The resulting activities for the individual 
chains and radio-nuclides, using the best fit of the data to the Monte Carlo simulations, are pre- 
sented in Table 0. The background of the detector above ~500 keV is thus currently dominated by 
the residual ^^^U, ^^^Th, '^''K and ^''Co activities in the Cu of the shield, and, to a lower extent, in 
the copper of the cryostat. Below ~500 keV, the background is dominated by ^'"Bi bremsstrahlung 
from ^^"Pb decays in the innermost lead shield. The obtained ^'°Pb activity is with 5.7ib0.5 Bq/kg 
higher than 3 Bq/kg, which is the value specified by the provider. 

5.2 Sample Analysis 

Prior to the screening of a specific sample, an estimate of the minimal measuring time tmin, based 
on the mass, shape and the targeted activity Atarget is performed: 



^target ■ r -m- £ 



(5.1) 



^ A background run up to an energy of ~8 MeV was acquired for a few weeks and revealed no U/Th contaminations 
close to the crystal, which would be visible as alpha peaks above a few MeV. 
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Table 5. Activities of residual and cosmogenic radionuclides inside the copper shield, the Ge crystal and the 
cryostat and of ^^^Rn inside the shield. 



Isotope 


Specifi 


c activity [^uBq/kg] 


Activity [^uBq/m^] 


Activity [Bq/kg] 


(chain) 


Cu (shield) 


Ge crystal 


Cu (cryostat) 


Sample chamber 


Inner Pb shield 


226Ra(238u) 


56 ± 11 




8±5 






228Th (232xh) 


27 ±7 




4±2 






40k 


32 ± 13 


<1.30 


11 ± 6 






60Co 


8 ±4 


<0.80 


1.3 ± 0.4 






''Co 


<0.27 


<0.11 


<0.22 








<1.30 


<1.60 


<2.15 








<0.16 


<0.15 


<0.50 






222Rn (238U) 








<55 




210pb 










5.7 ± 0.5 
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Figure 6. (Top) Comparison of the sum of the simulated spectra from natural and cosmogenic radionuclides 
in the detector and shield materials (black) with the observed background spectrum (red data points). (Bot- 
tom) The individual, best-fit contributions to the observed spectrum are shown: natural radioactivity in Cu 
(blue), cosmogenic radio-nuclides in Ge and Cu (green), ^^^Rn decays inside the shield (magenta) and ^lOpb 
decays in the Pb shield (yellow). 

where m is the mass of the sample, r is the branching ratio, e is the detection efficiency and the 
detection limit is defined in equation ([4.2|). Typical numbers for are 1-3 weeks in order to 
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achieve a sensitivity below a few mBq/kg. Figure |7| shows the sensitivity versus counting time for 
different sample masses using a specific gamma line (right) and for various radio-isotopes for a 
specific sample (left). 




5 10 15 20 25 30 5 10 15 20 25 30 

Counting time [days] Counting time [days] 



Figure 7. (Left) Sensitivity versus time for '^''K (black), the chain (using the ^''^Bi 609 keV line, red), 
the 232Th chain (using the ^'^Pb 239 keV line, green) and ''"Co (blue) for a 13.5 kg sample of PTFE. (Right) 
Sensitivity versus counting time for different PTFE sample masses: 4 kg (black), 8 kg (red), 12 kg (green) 
and 16kg (blue) using the 609 keV gamma line from ^'"^Bi. 

The activities of the screened materials are calculated using the analysis methods described 
in Section ^TT] and As evident from Table ^ where a few results from screened samples are 
provided, an excellent agreement between the results of the two methods was achieved. It provides 
a cross-check on the reliability of the used methods and on the obtained activities (or upper limits) 
for the samples that were screened. A large set of measurements is described and discussed in 
detail in [El. 



Table 6. Comparison of results obtained by the two different analysis methods using four different screening 
samples. 



Material 


Method 


Activity 




228Th 


40k 


60Co 


(amount, time) 














Copper 




mBq/kg 


<0.5 


<0.5 


<0.9 


0.24 ±0.08 


(18 kg, 20 d) 


y-lines 


mBq/kg 


<0.3 


<0.3 


<1.3 


0.24 ±0.06 


Stainless steel 




mBq/kg 


4.1±0.6 


1.4±0.4 


<4.6 


7.4±1.1 


(6.6kg, 6.8 d) 


/-lines 


mBq/kg 


4.3±0.9 


<1.8 


<5.7 


7.2±0.9 


Concrete 




Bq/kg 


17±2 


3.1±0.7 


53±4 


<0.6 


(35 g, 0.7 d) 


y-lines 


Bq/kg 


15±2 


3.8±0.8 


42±6 


<0.7 


Photomultipliers 


x' 


mBq/PMT 


<0.2 


0.20±0.09 


8±1 


0.6±0.1 


(22 pieces, 5.5 d) 


/-lines 


mBq/PMT 


<0.2 


0.18±0.06 


11±2 


0.6±0.1 



6. Summary 

The Gator screening facility at LNGS, described in detail in this paper, includes one of the world's 
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highest sensitivity HPGe spectrometers. It allows to measure large samples and to reach a sen- 
sitivity on specific activities below ~mBq/kg when using typical sample masses of a few kg and 
measuring times extending from one to several weeks. The integral counting rate of the detector 
in the 100-2700 keV energy region has decreased by more than a factor of 5 from (0.842ib0.005) 
counts/min at Soudan to (0.157 ±0.001) counts/min at LNGS. This decrease is mainly due to the 
improved shield and radon protection system, but also because of the decay of cosmogenic isotopes 
in the crystal and surrounding copper. The background of the facility, which was modeled using a 
detailed geometry of the spectrometer and its shield, and Monte Carlo simulations with Geant4, is 
dominated by U/Th/K decays in the copper of the cryostat and the shield, and from residual ^^''Pb 
in the innermost lead shield. This results from a comparison of background measurements with 
the simulations, by choosing the best fit of the MC spectra to the measured data. Two different 
data analysis methods show an excellent agreement when applied on a variety of samples. De- 
tailed screening results for the XENONIOO [Ql experiment are presented and discussed in [EOl . At 
present, the Gator facility is operated under stable conditions at LNGS. It is used to screen materi- 
als for the construction of the XENONIT experiment, for the next phase of the GERDA project, as 
well as for a future noble liquid dark matter seaixh facility, DARWIN [HPH . 
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